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1,0 INTRODUCTION AND RESEARCH GOALS 

The Large Area Crop, Inventory Experiment (LACIE) demonstrated the usefulness 
of Landsat data for crop identification measurements In support of area estimation 
without a ground sample stage, In particular* LACIE area estimation procedures 
were for estimation of small grain acreages In areas where ground data would not 
be available. The LACIE small grains procedures were end-of-season procedures. 

That Is, the LACIE procedures depended upon the full set of crop season acqui- 
sitions In order to perform the crop type labeling operation with sufficient 
accuracy. During the LACIE Transition Program, and In the subsequent AgRI STARS 
Project, there was a desire to not only extend the LACIE-Hke technology to crop 
types other than small grains, but also to more directly address the need for a 
continual through-the-season estimation capability, Through-the-season capability 
was most desirable to all potential users of the technology, since crop production 
Indications prior to the end-of-season would allow for more effective management 
and planning within the agrobusiness community and commodities market. 

The crop types selected for the extension of the LACIE-Ilke technology were 
corn and soybeans. The particular foreign areas of Interest for application of 
the technology were Argentina and Brazil. Therefore the work described In this 
report was directed toward corn and soybean crop types. The ultimate goal of the 
work was to develop corn and soybean crop type labeling procedures to be used In 
support of a through-the-season area estimation technology. It was envisioned 
that a through-the-season capability would be developed that would allow for the 
estimation of crop acreages at any time In the growing season. This may mean that 
more than one pre-end-of-season labeling procedure would need to be developed for 
use In deferent periods within the growing season, At the least, It was felt that 
the particular point In the season at which an estimate was desired would affect 
the operating parameters of any labeling procedure.. Therefore, a further goal of 
this research was to come to understand the dynamics of crop separability as the 
growing season progressed, 


2.0 GENERAL APPROACH 

There were three phases planned in the overall approach to developing through- 
the-season labeling procedures. These three phases were: 

First Phaser Determination and Description of Summer Crop Type 
Separation Characteristics Throughout the Growing 
Season 

Second Phase:! Development of Labeling Guidelines for Pre-End-of- 
Season Labeling of Corn and Soybeans 


it 
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Thlrd Phase: Development of Labeling Procedures to Support a Through- 
the-Season Area Estimation Procedure 


Briefly in the first phase— Determination and Description of Summer Crop Type 
Separation Characteristics Throughout the Growing Season— the research would focus 
on: determining what, If any, crop type or crop group characteristics could be 
observed within Landsat data that would be useful In uniquely Identifying the crop 
types or at the very least the crop groups of Interest. In addition, the expected 
Improvement In labeling accuracy as the growing season progressed was to be eval- 
uated and some mechanism for determining or for controlling the setting of the 
operating parameters i\t various periods In the growing season was to be Investi- 
gated. 

The second phase of research— Development of Labeling Guidelines— would form- 
ulate consistent labeling guidelines for the separation of corn and soybeans 
throughout the growing season. Clearly formulated guidelines would then allow 
labeling procedures to be designed In support of specific through-the-season area 
estimation procedures being developed. 

The third phase— Development of Labeling Procedures— was to procedural Ize the 
labeling guidelines Into a consistent logical processing of the data either by an 
analyst, machine, or both Interactively, similar to the procedural izatjon effort 
undertaken In the design of the Corn/Soybean Baseline Procedure C/S-l.' 

This document will report on the progress of the first phase of the research. 
At present no plans exist to continue with the other two phases due to termination 
of support for this effort. 


3.0 THROUGH-THE-SEASON MODEL OF LABELING ACCURACY 
3.1 GENERAL MODEL CONCEPT 

There are four characteristics observable in multi temporal Landsat data 
that are most useful in determining the crop type or crop group identification 
of an agricultural feature. These characteristics are: (1) a temporal char- 

acteristic which pertains to the development of a feature's vegetation canopy 
through time, i.e., what time of year is a green vegetation canopy observable 
for the feature as opposed to bare ground; (2) a spectral characteristic 
which refers to the value of a given spectral variable, or to the relative 
relationship between two or more spectral variables for a given feature; and 
(3) two spatial characteristics, one which refers to the geometric attributes 
of a given feature such as the nature of a feature's boundary, and the other 
which refers to the areal distribution of similar and dissimilar features 
within the area of analysis. The most important characteristic for crop 
identification is the temporal characteristic. The temporal characteristic 
allows the separation of crop groups within a region, or more precisely, is 
used to define the various crop groups within the region. In order to support 
evaluation of the temporal characteristic, collateral data in the form of crop 
calendars are needed. Indeed, the most Important piece of collateral data for 
crop type analysis is crop calendar data. The crop calendar data provides in- 
formation about when crops within a given region are likely to be or are (if a 
real time crop calendar Is used) in any given biostage. 


From past research and experience, the spectral response of many vegetation 
biostages are familiar and generally known. Therefore a spectral crop calendar 
can be generated from phenologlcal crop calendar data. This spectral crop 
calendar model allows an analyst or programmed decision algorithm to judge the 
most likely crop type or cover class label to assign to a given feature In a 
given set of mul tl temporal , spectral data such as Landsat data. 

The logic process by which a labeling target's observed characteristics are 
evaluated can be Illustrated by the list of "Important questions" In Table 1, 

The ability to answer the Important questions In a deterministic manner Is 
dependent upon the temporal set of acquisitions available for analysis. It can 
be appreciated that as the crop growing season starts, there Is little or no 
Spectral data available to answer any of the questions. As the season progresses, 
however, more data become available and more of the "Important questions" can be 
answered. The accuracy of any labeling analysis task Is a function of the tem- 
poral sequence of acquisitions available. Thus, It can be understood that the 
labeling accuracy Increases as a function of time. Note, however, that some 
questions can only be answered after a discrete event has occurred. For example, 
question 2a ("When does the maximum value for the green vegetation indicator 
occur?"): This question cannot be answered until some time AFTER the occurrence 

of the maximum value because a relative maximum can only be Identified after the 
value of the green vegetation Indicator has started to decrease again. This time 
delay in ability to answer a specific question relative to the occurrence of some 
events says that the labeling accuracy function relative to time Is not a 
smoothly increasing function, but is a function which increases In steps. The 
plot of this function through the crop growing season would have several plateaus 
where the labeling accuracy would remain fairly constant until sufficient data 
had been received that would allow the answering of additional "Important 
questions," (A graphical representation of this model is shown In Figure 1.) The 
above description of the model of labeling accuracy has ignored the question of 
data dropout and the use of. multiyear labeling procedures or estimators. If multi 
year labeling procedures or estimators were employed, the value for A Q (initial 
value for labeling accuracy at the start (t 0 ) of a given crop year) would probably 
be higher (indicated by A,,,). The labeling accuracy function would then proceed 
to perform generally as described above except for the increased accuracy levels 
expected from the use of multiyear data. Data dropout would have the effect of 
prolonging the length of various labeling accuracy plateaus until such time as 
sufficient data became available to answer additional "important questions." 

With the above-described Labeling Accuracy Model as a basis, the Through-the- 
Season label ing research effort sought to accomplish the following objectives in 
Phase One: 

(1) Determine the specific identification characteristics that could be 
used to uniquely identify summer crop types (corn/soybean) or crop 
group to the highest accuracy level possible within any given accuracy 
plateau; 

(2) Determine a means of specifically identifying when each of the accuracy 
plateaus started and ended for the corn/soybean situation. 


3.2 SPECIFIC THROUGH-THE- SEASON MODEL OF LABELING ACCURACY FOR 
CORN/SOYBEANS IN U.S. CORN BELT AND PERIPHERY 

The accuracy plateaus described In the general Through-the-Season labeling 
model need to be more definitively described for the specific corn/soybean 
situation. This was done based on when each of the Important questions from 
Table 1 could be answered during the Corn/Soybean Summer Crop growing Season 
In the U.S. Corn Belt and Corn Belt Periphery. Figure 2 shows the specific 
characterization of the Labeling Accuracy Model for the Corn/Soybean Summer 
Crops In the Corn Belt Area. Use of multiyear data or data dropout were 
Ignored for the present, The C/S-l (Corn/Soybean Baseline Procedure) Summer 
Crop biowindows and the potential and actual separation windows have also 
been delineated on the specific model representation. As can be seen, C/S-l 
Is! actually a mid-season procedure based on the occurrence of the actual 
separation window In the mid-season of the model. Since the DFS strat- 
ification In C/S-l can be done with a minimum set of two acquisitions— one 
acquisition from the very early or pre-season time period and one acquisition 
from the very early to mid-season time per1od--C/S-l can operate as a mid- 
season procedure. Partly because of the above reason and partly due to 
efficiency of approach, the decision was made to concentrate on the early 
season plateau In the Initial efforts toward developing a Through-the-Season 
labeling logic. 

Therefore the specific objectives of this Initial Through-the-Season 
effort were: 

(1) to determine when t 2 occurred - for the Corn/Soybean situation, 

(1) Investigate how the occurrence of t£ could be predicted 
operationally (l.e. real time during an actual Through- 
the-Season Estimation Effort) 

(2) to determine what characteristics or attributes of corn/soybeans, 

summer crops, and others were useful In separating 

(I) summer crops from other (at the least) 

(II) corn and soybeans from other summer crops and from 
each other if possible, 


0 DATA ANALYSIS METHODS 
4.1 DATA SET 

4,1.1 Segment Selection 

The Through-the-Season labeling guidelines research was intended 
to support Through-the-Season area estimation procedures for use in 
foreign areas, particularly Argentina. However at the time of this 
study, no Argentina spectral data were available. For this reason 
the labeling guidelines research was conducted using the spectral 
data from segments located within the U.S. Corn Belt and Corn Belt 
Periphery. Segments selected for Inclusion in the study were chosen 


based on the following; (1) how well they represented the geo- 
graphical diversity within the Corn Belt and Corn Belt Periphery, 
and (2) the adequacy of their Landsat acquisition histories, A 
good acquisition history was one that Included at least one ac- 
quisition prior to the earliest planting date for corn, at least 
one acquisition during or just past corn tassellng or soybean pod 
setting, and a sufficient number of acquisitions In between these 
two bounding crop biostages. The dates of corn planting, corn 
tassellng, and soybean pod setting, as well as other crop calendar 
Information was determined from the statewide year-specific crop 
calendars published by the U.S.D.A. - Statistical Reporting Service 
within their Weekly Crop and Weather Bulletins. Before final In- 
clusion of a segment within the study data base, all segment ac- 
quisitions were visually screened for cloud cover, haze, and bad 
data lines. Eighteen segments were selected to comprise the data 
set. The list of selected segments and acquisitions Is given In 
Table 2, and the segment locations are shown In Figure 3. 

4.1.2 Segment Preprocessing 

Sun angle and haze correction was performed on all acquisitions 
selected for Inclusion in the data set. The sun angle and haze 
correction algorithm used to normalize the spectral data was the 
global XSTAR algorithm 2 developed by the Environmental Research 
Institute of Michigan (ERIM). This algorithm had been used quite 
Successfully on numerous previously analyzed data sets of Landsat 
spectral data. Following the application of XSTAR, the spectral 
features of Interest for the analysis were extracted. The spectral 
features extracted for analysis were as follows; 


(a) Tassel-cap component 1 (TC-1, Brightness) 

■(b) Tassel-cap component 2 (TC-2, Greenness).- -Green 
vegetation Indicator 

(c) GRABS 3 (Greenness Above Bare Soil)— Green vegetation 
indicator 


Tassel-cap Brightness (TC-1) and GRABS were given the most emphasis 
in the analysis due to the successful results achieved in labeling 
crop type from prior studies. 

Tassel-cap transformation is an afine transformation of the landsat 
four-channel MSS data into four new channels or components of data 
which correspond more directly with physical processes and relation- 
ships, (particularly for vegetational ground features) than the 
Landsat MSS bands do directly. GRABS is a green vegetation Indicator 
(measure of infraredness) derived from tassel-cap greenness and the 
threshold of detection for green vegetation established for the 7/5VI 

GRABS * (TC-2) - 0.09178 (TC-1) + 5.58959. 


4.1 .3 Sampling of Fields by Crop Type 

Digitized ground data supplied by Johnson Space Center (JSC) 
was used to Identify specific fields within a given crop type or 
vegetation cover class. A maximum of 10 fields for each summer 
crop type present within a segment was sampled. Sufficiently 
large non-summer crop fields were also sampled If they were 
present within a given segment. To minimize the effects of mis- 
registration between acquisitions, and to Insure some level of 
spectral purity within the crop type samples, only field center 
pixels were Included within the sample. A pixel was Included In 
the sample of a field If It had a margin around It that was at 
least two pixels wide, or If the surrounding pixels were of the 
same cover class. This sampling strategy had the effect of ex- 
cluding small fields from the study. In total 574 fields were 
sampled from the 18 segments In the data set. The number of 
fields sampled by segments and cover class Is shown In Table 3. 

4.1.4 Analysis Approach 

For each of the sampled fields, a field mean was calculated 
In each of the spectral features to be analyzed for each acqui- 
sition. All of the analyses were then carried out on the field 
means data. Temporal plots of the field means were generated for 
Brightness (TC-1) and GRABS spectral features for each sampled 
field. In addition, segment level (segment average) temporal 
plots were generated for each segment by using a weighted (by 
pixel number) average of the field means within each crop type 
within each segment. Appendix A contains the segment mean tem- 
poral plots for Brightness and GRABS for each cover type sampled, 

In order to determine whether there were any detectable 
differences between the crop groups, and the summer crop types 
at various points earlier in the growing season prior to end-of- 
season, the above described temporal plots were examined and eval- 
uated with particular emphasis on observing differences as early 
In the growing season as possible. In addition the data were eval- 
uated to try to determine at what points In the growing season 
significant increases in identification accuracies would likely 
occur as a function of the additional spectral and crop calendar 
data which become available as the season, progressed. 


5.0 DETERMINATION OF t 2 : SUMMER CROP SPECTRAL EMERGENCE DATE 

5,1 DEFINITION, RATIONALE AND ASSUMPTIONS 

In the Corn Belt, soybean is frequently the last summer crop planted. 
This is so partly because corn yield is more favorably affected by early 
planting than is soybeans. Thus corn planting has priority over soybeans 
also soybeans can be double cropped after small grains so that double 
cropped soybean planting Is delayed until after small grain harvest. 
Observations from past research of the temporal patterns of corn, soybean 


7- 


and other summer crops In the Corn Belt have shown this later planting 
pattern for soybeans to be reasonably consistent within the region, 
though local variations do occur. The date of tg Is defined as the 
date when all summer crops are detectable. Within the Corn Belt tg can 
be approximately defined as the date when all soybeans are detectable. 

The word 1r all" In these definitions should be read >"15%", This Is 
similar to the manner In which the spectral biowindows were defined for 
C/S-l, 

Since by definition from the model, the ’'early season" labeling 
guidelines would be based on the spectral -temporal characteristics of 
detectable vegetation, It was necessary to know or to be able to est- 
imate when "all" summer crops were detectable In order to correctly 
operate the appropriate labeling procedures. * That Is, the actual date 
of occurrence of tg needed to be determined, A method for estimating 
the occurrence of t2 prior to its occurrence was also expected to be 
needed In the eventual Thrpugh-the-Season Estimation Procedure to be 
developed. 

Since t 2 was expected to vary between years, a long-term average 
crop calendar would not be adequate for estimating t 2 as closely as 
was felt necessary. Therefore, a real-time method for estimating the 
occurrence of t 2 In any given Inventory year was needed, The assumption 
was made that, at the least, some real-time crop calendar model data would 
be available. However, since no crop calendar model data was available 
for evaluation, and only limited progress on a Corn Crop Calendar Model 
was perceived at the time by these Investigators, reliance upon such crop 
calendar model data was minimized as much as possible. The only year- 
specific crop calendar data available for this study, however, was the 
U.S.D.A.-S.R.S, 's Statewide Crop Calendars, Therefore the statewide year- 
specific crop calendar Information was used to develop the method for est- 
imating the occurrence of t 2 , It was hoped that crop calendar models under 
development could supply the needed Information In foreign areas at some 
time In the not-too-dlstant future. 

Since t 2 was expected to vary each year within each region, a ref- 
erence point or start date tg was required that could be determined for 
each region each year. From the statewide year-specific U.S.D.A.-S.R.S, 
Crop Calendars tg was defined as the day on which 50% of the first planted 
summer crop was reported to have been planted.. Thus In the Corn Belt, tg 
was defined as the date on which 50% of corn was reported to have been 
planted. From past research In crop calendar analysis for summer crops 
and small grains, It was determined that the 50% planted date for any crop 
tended to be less variable year-to-year than other percent planted points 
e.g. 5%, 75%, 95%, etc. Based on this perceived tendency toward higher 
stability, the 50% planted date was chosen as tg. In addition, this def- 
inition of tg allowed tg to be determined from data that would be avail- 
able reasonably early In any growing season so that true Through-the-Season 
procedures could be developed that would start with the season and not at 
some later point in the season. (Past definitions of early and mid-season 


Labeling guidelines for a period <t 2 are called Very-Early Season 
guidelines by definition and may not rely on detectable vegetation 
for all potential summer crop fields, 


h»vt been stated In terms of days preceding harvest* One, therefore, 
needed an estimate of when harvest was to occur before one could 
classify acquisitions as to early or mid-season.) 

The time between t 0 and t 2 (time Interval t 2 -t 0 ) needed to be 
determined and Its year-to-year variation evaluated. 

5.2 ESTIMATION OF (t 2 -t 0 ) 

Figure 4 graphically represents one preliminary method used to 
estimate (t 2 -t 0 ). 1978 data from segment 127, Montgomery County, 

Indiana, was used to make a "first-look" sample calculation of the 
Interval from 50% corn planted (to) to 95% soybean detectaulllty (tz). 
Segment 127 was selected from the 18 sampled segments because It had a 
good acquisition history, good crop calendar Information, and was one 
of the first segments to pass through pre-processing (sun angle and 
haze correction, etc.). GRABS values were used to calculate the detect- 
ability u'ate because GRABS»0 (green vegetation threshold value) Is 
standard for all segments and acquisitions processed through She LACIE 
Program. Mean GRABS values of ten soybean fields were plotted for the 
first five acquisitions; June 1, June 10, July 16, July 26 and August 4, 
All ten fields passed fh si GRABS threshold of detectability between the 
fourth and fifth acquisitions, Figure 5. A variable x (Julian date) « 
the t Intercept, of the line drawn between GRABS values on acquisitions 
4 and 5 with the threshold of detectability (GRABS * 0), was determined 
for each sampled field. Using the value of this variable for each sample 
field, a weighted mean and variance was calculated. The 95% confidence 
Interval was determined and added to the weighted mean ving a Julian 
day by which 95% of sampled soybean fields were detecUule. In 1978 50% 
statewide planted day for soybean was obtained from the Indiana annual 
crop and livestock report for the year and subtracted from the 95% mean 
detectability day + the 95% confidence interval. The difference, 15 days, 
was the Interval x 2 In Figure 4, 

The Interval from the 50%-pl anted day for corn and the 50%-pl anted 
day for soybean was ascertained as shown In Table 4. Six years of crop 
calendar data was available from which to calculate a mean number of days 
from the statewide 50% corn planted day to the 50% soybean-planted day. 

The 95% confidence interval was determined and added to the mean to give 
an interval, xj in Figure 4, of 26 days. The 26-day Interval was added 
to the 15-day interval from soybean planting to first detectability to 
estimate x^j in Figure 4— the number of days from the 50%-pl anted day for 
corn to first detectability of all soybean fields in segment 127. The 
final value, 41 days, was rounded off to 40 days for convenience. 

Although the 40-day Interval from corn planting to soybean detect- 
ability was calculated from a single segment year, the time interval 
tg+40 was used as the initial estimate of (t 2 -t 0 ) for all 18 segments. 

Tne actual value for t 2 -tg varied from state to state because of state- 
to-state variation in the 50% corn-planted day, but almost all corn and 
soybean fields were detectable on the first acquisition after tg+40. 


99.1% of sniper crop fields had mean GRABS values above 0 on the first 
acquisition after tg. 97.3% of the fields had mean "slues above 
GRABS ■ 5.0* The number of days from tg+40 to the first acquisition 
varied between 5 and 40 with an average, weighted by the number of 
fields In each segment* of 18.5 days. In contrast, 41.8% of summer 
crop fields had GRABS values above 0 on the last acquisition prior to 
the tg+40 date. The number of days between the acquisition prior to 
tg+40 and the (tg+40) date varied between 5 and 31 with a weighted 
average of 15.0, The data Is summarized In Table 5, 

The estimate of a 40-day Interval from planting to detectability 
was further evaluated using field measurements data obtained from LARS, 
Corn fields had at least 20% ground cover* approximately 40 days after 
planting, especially with the more dense, and more common, plantings. 
Soybean fields passed the 20% ground cover point well before 40 days 
after 50% planted date of the soybeans. The spectral emergence of the 
soybean fields could not be evaluated with respect to a 50% corn 
planted date since this was somewhat of an artificial situation. 

Additional analysis was performed to estimate t 2 using all corn and 
soybean pixels In seven Indiana and seven Iowa segments. Digitized and 
registered ground data was used to sample all Identified corn and soy- 
bean pixels within the 14 segments. The percentage’ of corn and soybean 
pixels above GRABS ■ 5.0 decision threshold** on each acquisition for 
each segment was determined and plotted against Julian date. The seg- 
ment plots were then combined by state and a statewide (as represented 
by the sample of seven segments/state) cumulative spectral emergence 
curve (S curve) was generated for each state. The S curve was fit 
to the data points using a two piece curve fitting method whereby all 
data before the obvious upper break In slope plus the 1st group of data 
points at or after the obvious break In slope were fitted by* regression 
to the exponential function y * ae bx . The second group of points which 
Included all data points at or after the obvious upper break In slope 
were fit by regression to the linear function y » ax+b. The two 
curves were then pieced together visually (i.e, the upper shoulder was 
estimated by eye). Two analyses were done, a) one using tg * 50% planted 
date for corn statewide, and b) one using tg ■ 50% planted date for corn 
at the CRD to which a given segment belonged. Thus the second analysis 
was done on data normalized for tg at the CRD level, The non-normal 1 zed 
S curves of data points for each set of segments by state are shown in 
Figure 6. The normalized S curves and data points are shown In Figure 7, 

The results of this analysis are shown In Table 6. The earlier est- 
imate of tg+40 * t 2 ■ 95% summer crops spectral emerged based on analysis 
of segment 127 sampled corn and soybean fields does not appear to have 
been very good. From the above described analysis t 2 * tg+70 and tg+51 
for Iowa and Indiana non-normal 1 zed data respectively. For the normalized 


From work at ERIM it has been estimated that approximately 20% ground 
cover Is necessary In order to reach the threshold of spectral detect- 
ability. 

** Two GRABS detectability thresholds were used, 0 and 5, because temporal 
plots of summer crop field mean values indicated that GRABS values some- 
times fluctuated between 0 and 5 prior to tg+40. GRABS values did not 
exceed 5 on any consecutive acquisitions prior to tg+40. 
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data , t 2 ■ To+58 and Tq+ 56 days for Iowa and Indiana respectively. There 
was also better correspondence between estimates of other t emergence 
points between the two states In the normalized data. The r 2 values for 
the regression curve fit were fairly good for Indiana - .86 and 44 fgr 
non-normal 1 zed and normalized data respectively. However the r 2 yatMS 
for the regression curve fit for Iowa were not very good - #57 and 
non-normal 1 zed and normalized data respectively. 

This new estimate for tg 8 to+57 days actually corresponds very well 
with the definition for the end of Summer Crop biowindow SC2 from C/S-l 
baseline procedure If one adjusts for the difference between 5 % corn 
planted (from C/S-l) and 50* corn planted from the Through-the-Season 
Model. This was to be expected since the SC2/SC3 biowindow boundary 
point and t 2 are theoretically the same from their conceptual definitions. 

It Is clear from Inspection of Figures 6 and 7 that there is much 
more variation In Iowa spectral crop calendars. To what this variation 
Is due cannot be determined In any certainty In this study, however It 
wOuld seem reasonable to suspect that most of the crop calendar variation 
In Iowa Is due to greater variation In planting dates throughout the state 
than occurs In Indiana. The observation that the two normalized fitted 
spectral emergence curves for Iowa and Indiana were similar leads one to 
suspect that normalization for planting date variation can be valuable for 
any agricultural multitemporal spectral data processing. This of course 
was the Idea behind UCB's dynamic stratification work. However the large 
variation In the Iowa data even after normalization seems to suggest that 
crop calendar "strata" In Iowa are smaller than CRD size, In Indiana on 
the other hand CRD size strata seem to be reasonable. 

The above observations are relevant to the Through-the-Season effort 
In that It suggests that any Through-the-Season Labeling Procedure or 
Area Estimation Procedure based on certain kinds of models (l.e. models 
that have some variable based on crop calendar) would be more effective 
If normalization for planting date can first be performed. Another method 
for controlling the variation would be through crop calendar stratification. 

Which method would be more desirable would be totally dependent upon 
the over-all area estimation procedure design and the operational environ- 
ment of the procedure. 


6.0 EARLY SEASON TEMPORAL-SPECTRAL CHARACTERISTICS: SUMMER CROPS VS OTHER 

Edwin J. Sheffner 

6.1 GREEN VEGETATION INDICATOR-GRABS 

Differences between summer crops and non-summer crop temporal profiles 
become more apparent when examined within the context of the Through-the 
Season Labeling Accuracy Model where t 2 is defined as a reference point. 
Note: The following analysis was performed prior to the latest estimation 
of t 2 - to+57 days, and so the analysis was performed relative to the 
original estimate of t 2 ■ tg+40 days. 
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6.1 .1 Summer Crops vs. Other 
Spectral Emergence 

Mean GRABS values for corn, soybean and sorghum fields rarely 
exceeded the decision threshold GRABS ■ 5.0 on acquisitions prior to 
the tg+40 date. 1572 was the only segment with summer crop mean 
GRABS values above 5.0 on any acquisition prior to tg+40 (sorghum), 
and no segments had two consecutive acquisitions prior to tg+40 with 
summer crop values above 5.0. Sixteen of the 18 sampled segments 
had alfalfa /pasture fields. Mean GRABS values of alfalfa/pasture 
exceeded 5.0 on at least one acquisition prior to Y Q +40 In all 16 
segments. Fourteen of the 16 segments with alfalfa/pasture fields 
had more than two acquisitions prior to tg+40; In 12 of the segments 
alfalfa/pasture field mean GRABS values exceeded 5.0 on at least 
three acquisitions prior to tg+40. 

Idle cover and fallow fields also tended to have higher GRABS 
values than summer crops early in the season. Eight segments had 
sampled Idle cover fields and six segments had sampled fallow fields. 
Idle cover field mean GRABS values exceeded 5.0 on more than one 
acquisition prior to tg+40 In five of the eight segments and exceeded 
5.0 on more than two acquisitions prior to tg+40 in four of the eight 
segments. Fallow field mean GRABS values exceeded 5.0 on more than 
two acquisitions prior to tg+40 in three of six segments. 

The sampled summer crop confusers tended to reach maximum GRABS 
values earlier in the season than summer crops. All summer crop 
fields reached maximum GRABS values after the tg+40 date, but alfalfa/ 
pasture fields reached peak GRABS values prior to tg+40 In 12 of 16 
segments. Idle cover fields reached peak GRABS values prior to tg+40 . 

In two of six segments. 

A dip in summer crop mean GRABS values prior to the tg+40 date 
occurred in several segments. The GRABS temporal plots or corn and 
Soybean In segment 127, (Figure 5), illustrate the dip in GRABS -values 
which occurred between the third and fifth acquisitions. The CRABS 
value dip was a more common occurrence in soybean and sorghum fields 
than in corn fields. The dip occurred in soybean fields in 14 of 17 
segments and in three of four segments with sorghum, but only in half 
the segments (9 of 18) with corn. Similar dip phenomenon was noted in 
non-summer crops. Nine of 16 segments’ with alfalfa/pasture, two of 
eight segments with idle cover fields and two of six segments with 
fallow fields had similar dips prior to tg+40. The cause of the dip 
is Uncertain. It may be a response to cultivation practices or weather 
or variation in sensor calibration. In either case, Its appearance in 
the spectral data is acquisition dependent. 

6.1.2 Differences Between Summer Crops 

The GRABS temporal plots showed differences among the summer crops 
In timing and range of peak GRABS values, Corn fields generally peaked 
earlier than soybean fields and sorghum fields were mixed. Mean corn 
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field values reached a peak between t 0 +4Q and tn+85 In 11 of 17 
segments , whereas soybean fields reached peak GRABS values after 
date tn+85 In 10 of 15 segments. Three segments lacked adequate 
acquisitions to use In this analysis. Segments with sorghum fields 
were evenly split— two of the four segments had sorghum fields which 
peaked In the period t 0 +40 and tn+85 and two segments had sorghum 
fields which peaked In the period after to+85. During specified 
periods of the growing season, mean summer crop GRABS differed. In 
the period between tg+40 and to+57, corn mean GRABS values varied 
between 14.42 and 31.12. Mean soybean field values were lower, 
varying between 1.03 and 27.62. Later In the season, the situation 
reversed. For acquisitions obtained between tn+85 and to+102, mean 
corn values were between 17.62 and 22.63. Soybean mean field values 
were between 31.52 and 45,65 during the same period. Sorghum field 
GRABS values varied between 7.17 and 12.18 In the earlier period and 
between 15.67 and 31.12 later In the season. 

GRABS emerged as a significant feature for early season separ- 
ation of summer and non-summer crops. Vegetation within fields of 
major summer crop confusers tend to be detectable on Landsat In the 
GRABS feature earlier than summer crops, reach peak values before 
summer crops, and generally described different spectral development 
patterns, The descriptions of the development of the various cover 
types given above and the conclusions drawn from them are based on 
mean values of all field means of a given cover type. Individual 
fields varied significantly from the mean as can be seen In 
Appendix A by the plots of standard deviation. Some non-summer crop 
fields, particularly Idle cover and fallow fields* mimic sunnier 
crops In GRABS development pattern, and, occasionally, a field ident- 
ified as summer crop on ground data was observed to follow a typical 
alfalfa/pasture pattern, The possibility that anomalies are present 
must be kept In mind during the labeling process. 

6.2 BRIGHTNESS VALUES 

The dip noted In the GRABS feature for summer crops Is also observed 
In the tassel-cap brightness feature. Mean TCI values for corn fields In 
12 of 17 segments dipped and then Increased In value prior to tg+40. The 
same behavior occurred for soybean fields in 14 of 16 segments and 3 of 4 
segments with sorghum. The TCI dip also occurred in some non-summer crop 
fields, Mean values for Idle cover fields In five of eight segments and 
for fallow fields In 4 of 6 segments had a> TCI dip, but the dip was evident 
In only 5 of 15 segments with alfalfa/pasture. 

The time of maximum TCI values varied among the summer crops and be- 
tween the summer crops and non-summer crops. Maximum TCI values occurred 
prior to t 0 +40 In corn fields in 13 of 17 segments and In 3 of 4 segments 
with sorghum fields, but only In 4 of 16 segments with soybean. Maximum 
soybean TCI values were attained between t 0 +40 and to+85 In 6 segments and 
after t 0 +85 In 6 additional segments. 
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Maximum TCI values were reached In half the segments with alfalfa/ 
pasture fields prior to t 0 +40 and between t o +40 and t 0 +85 In the 
majority of the remaining alfalfa/pasture segments, 

TCI values for soybeans remained consistently higher than corn 
or sorghum values after the tn+40 date. Soybean field mean values of 
TCI varied between 56.38 and 81.68 on acquisitions between t 0 +40 and 
tn+57 and between 75.10 and 90.74 for acquisitions between tg+85 and 
t 0 +102. Corn field values for the same periods varied between 61 .73 
and 74.98, and 6179 and 83.76 respectively, Sorghum TCI values were 
similar to corn In both time frames. 


7.0 SUMMARY 

Landsat data has proven to be effective In large scale agricultural Inven- 
tories. However, prior to AgRISTARS, the usefulness of Landsat data was limited 
by the lack of labeling logic for Identification of crops other than small grains 
and earlier than at-harvest, 

To Improve the timeliness of Landsat area estimates, UCB under took the 
experiment described In this report. The results of the experiment Included a 
framework for understanding how Landsat data can be used In conjunction with 
ancillary data to Identify crop groups (In the Corn Belt) at any time during the 
growing season (labeling accuracy plateaus) and the demonstration that first 
detectability of green vegetation on Landsat was a critical event In the spectral 
development of summer crops in the Corn Belt which could be used to separate 
summer crops from non-summer crops as early as 40 days after the 50^-planted date 
statewide for corn (tQ+40), Maximum green development was Identified as the 
earliest time in the growing season when corn and soybeans are separable. 

While the results of this experiment indicate that procedures for summer crop 
labeling prior to harvest are feasible, there are a number of significant problems 
that remain to be examined. The method described in this report for determining 
the boundaries of the labeling accuracy plateaus assume that a season-specific 
starting date (from crop calendar information or a planting model) will be avail- 
able, The assumption may not be realistic. 

Although to+40 may have no meaning outside the Corn Belt, the labeling accuracy 
plateau appraoch may be appropriate In other regions. The nature of summer crop 
cultivation in non-Corn Belt regions needs to be studied before Through-the-Season 
labeling procedures can be written. 


Within the Corn Belt, additional research is required to quantitatively est- 
ablish the remaining boundaries of the labeling accuracy plateaus. A multiyear 
approach which Incorporates Information on crop rotation and planting patterns may 
yield a reliable date prior to to+40 for summer crop/non-summer crop discrimination. 
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TABLE 1. List of "Important Questions" for Crop Type Identification 



Question Statement 

Determinant 

Characteristic 

1 

When (during what time period) Is a green 
vegetation canopy detectable or present within 
the given feature? 

Temporal 

OR THREE MORE SIMPLE, ALTERNATIVE QUESTIONS: 


la 

Is there a bare soil stage observable? When? 

Temporal 

1b 

When Is a green vegetation canopy first detect- 
able within the spectral data for a given 
feature? 

lc 

When Is a green vegetation canopy last detectable 
within the spectral data for a given feature? 

Temporal 

2 

When do the maximum (or minimum) values for 
specific spectral variables occur? 

Temporal 

2a 

When does the maximum value for the measure 
of "infraredness the green vegetation 
Indicator, occur? 

Temporal 

2b 

When does the maximum or minimum value for the 
brightness measure (Brightness) occur? 

Temporal 

2c 

When does the maximum or minimum value for the 
visual yellow-green measure occur? (Can not 
be adequately measured with Landsat at this 
time) 

Temporal 

3 

What is the value of the maximum for the Green 
Vegetation Indicator? 

Spectral 

4 

Is there a senescent biostage observable for 
the feature? 

Spectral 

5 

If there is a senescent biostage observable 
for the feature, when does it occur? 

Temporal 

6 

What is the shape of the temporal 'profile 
(Temporal plot) of the green vegetation 
indicator vs. time? (Does the profile have 
a plateau, is it broad-peaked, narrow- 
peaked, tall, squat, etc.?) 

Temporal -Spectral 
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TABLE 2, Location and Acquisition Dates - Data Set Segments 


Segment Location Acquisition Dates 


123 

Hamilton Co. 
Indiana 

89 

107 

152 

161 

197 

233 

269 



127 

Montgomery Co. 

89 

107 

152 

161 

197 

207 

216 

243 

252 


Indiana 

269 









141 

Madison Co. 
Iowa 

86 

103 

130 

166 

167 

212 

220 



144 

Wapello Co. 
Iowa 

130 

165 

183 

219 

220 

238 

246 

264 

274 

205 

Clark Co. 

93 

101 

137 

155 

164 

209 

218 

246 

272 


Missouri 

282 

290 








809 

Ogle Co. 
Illinois 

101 

163 

164 

209 

218 

244 

254 



832 

Adams Co. 
Indiana 

88 

97 

151 

160 

178 

232 

268 



843 

Henry Co. 
Indiana 

88 

97 

151 

152 

160 

178 

197 

251 


852 

Randolf Co. 
Indiana 

88 

97 

151 

160 

178 

232 

250 

268 


853 

Randolf Co. 
Indiana 

88 

97 

151 

160 

178 

232 

250 

268 


860 

Wells Co. 

88 

97 

107 

116 

151 

152 

160 

161 

178 


Indiana 

197 

232 

251 







864 

Crawford Co. 

87 

96 

141 

150 

159 

186 

222 

231 

248 


Iowa 

258 









865 

Crawford Co, 

87 

96 

141 

153 

159 

168 

186 

231 

249 


Iowa 

267 









880 

Monona Co. 
Iowa 

87 

96 

141 

150 

186 

204 

222 



881 

Monona Co. 

87 

96 

141 

159 

186 

213 

222 

231 

249 


Iowa 

267 









882 

Palo Alto Co, 

86 

96 

131 

141 

150 

159 

186 

213 

221 


Iowa 

231 

258 

267 

293 






1572 

Custer Co, 
Nebraska 

134 

153 

170 

189 

206 

224 

243 



1591 

Webster Co. 

134 

152 

169 

188 

205 

223 

241 

259 



Nebraska 
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TABLE 3: Number of Fields Sampled From Data Set 



1 

i 


i;} 


TOTAL: Summer Crops * 338 Non-Summer Crops a 236 

C * Corn A * Alfalfa 

SY 3 Soybean P * Pasture 

SR 3 Sorghum H = Hay 

SB * Sugar Beet IC * Idle Cover 

SU 3 Sunflower F ■ Fallow 

T * Trees 
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Table 4: Calculation of maximum number of days between 50% planted 
day of Corn and 50% planted day of soybean In Indiana. 


Year 

A 

Day on which 
50% of Corn 
was Planted* 

B 

Day on which 
50% of Soy 
was Planted* 

B-A 

1973 

140 

161 

21 days 

1975 

133 

140 

7 

1976 

124 

140 

16 

1977 

131 

142 

11 

1978 

152 

160 

8 

1979 

138 

148 

10 


136.33 +24.48 

148.5 + 25.07 



Mean a 12.17 days 
Standard Deviation ■ 5.34 days 
NOTE: 95% Confidence Interval * 12.17 days + 13.74 
Day Interval » 12.17 + 13-.74 * 25.91 or 26 days 

* From Indiana Annual 
Crop and Livestock 
Summaries 
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TABLE S: Accuracy to t 0 + 40 Date In Predicting Summer 
Crop Detectability on Landsat 


Segment 

t 0 + 40 Day* 

A 

B 


C 


D 

123 

192 

5 

12 

12 

(100) 

12 ( 

|ioo) 

127 

192 

5 

20 

20 

(100) 

20 

!ioo) 

141 

174 

38 

13 

13 

(loo) 

13 

[100) 

144 

174 

9 

15 

15 

(100) 

15 

[100) 

205 

178 

31 

21 

21 

(100) 

21 

[100) 

809 

184 

25 

20 

20 

(100) 

20 

:ioo) 

832 

192 

40 

20 

20 

(100) 

20 

[100) 

843 

192 

5 

20 

20 

100) 

19 

[95.0) 

852 

192 

40 

20 

20 

(100) 

20 

!ioo) 

853 

192 

40 

20 

20 

(100) 

20 

[100) 

860 

192 

5 

16 

16 

(100) 

15 

[93.8) 

864 

174 

12 

16 

16 

(100) 

16 

(100) 

865 

174 

12 

24 

24 

(100) 

24 

[100 

880 

174 

12 

23 

23 

(100) 

23 

[100) 

881 

174 

12 

20 

19 

(95.0) 

18 

(90.0) 

882 

174 

12 

20 

20 

(100) 

20 

(100) 

1572 

175 

14 

10 

8 

(80.0) 

8 

(80.0) 

1591 

175 

13 

20 

20 

(100) 

20 i 

(100) 


A * Number of days from t 0 + 40 date to first acquisition after t 0 +40. 

B ■ Number of sampled summer crop fields. 

c * Number of sampled summer crop fields with. mean GRABS values above 

GRABS =0 on first acquisition after t 0 + 40 (%). 

D * Number of sampled summer crop fields with mean GRABS values above 
GRABS*5 on first acquisition after t 0 + 40 (%), 

* * t 0 + 40 Day * 40 days after 5056 corn planted day, statewide. 

Summary of Above Data 
Number of Segments 


4 

5 

68 

68 

(100) 

66 

(37.1) 

1 

9 

15 

15 

(100) 

13 

(86.7) 

5 

12 

103 

102 

(99.0) 

101 

(98.1) 

2 

13/14 

30 

28 

(93.3) 

27 

(90.0) 

6 

20+ 

114 

114 

(100) 

114 

(100) 
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TABLE 6: Evaluation of (t 2 - t 0 ) using accumulated 's-curve* approach. 



l * emerged 

- .todays) 

** emerged - 

todays) 

* Corn + Soybean 
Spectrally Emerged 

Non-Normal 1 zed 
t 0 ■ statewide 50* 
planted date 
of corn 

Normal 1 zed 
t fl ■ CRD 50* planted 
date of corn 

95* 

IOWA 
70 days 

INDIANA 
51 days 

IOWA 
58 days 

INDIANA 
55 days 

90* 

57 days 

*45 days 

52 days 

52 days 

75* 

49) days 

41 days 

47 days 

48 days 

50* 

42 days 

33 days 

39 days 

40 days 

25* 

29 days 

17 days 

25 days 

25 days 


Coefficients from y ■ ae bx curve fit 


a 

6.95 

10.38 

7.22 

7.91 

b 

.05 

.05 

.05 

.05 

r 2 

.57 

.86 

.46 

.84 
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ORIGINAL RAGE » 

OF POOR QUALITY 

FIGURE 3 

LOCATION OF SAMPLE SEGMENTS 



Figure 4 

Method for Determining when Summer Crops Have Passed 
Landsat Threshold of Green Vegetation Detectability 
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APPENDIX A 


✓ 

GRABS and Brightness versus Time Plots 


Plots of GRABS (.A-l to A-36) and Brightness tA-37 to A-72) 
versus time for corn, soybeans and other crops for sample seg 
ments listed In Table 2. The actual values equal the plotted 
value minus 10. The horizontal line Indicates the detection 
threshold (GRABS-0). The first vertical line equals t + 40; 
the second vertical line equals t + 85. 
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